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ABSTRACT 

The equations of motion for structures with adaptive elements for vibration control 
are presented for parallel computations to be used as a software package for real-time 
control of flexible space structures. A brief introduction of the state-of-tfc vart parallel 
computational capability is also presented. Time marching strategies are developed for 
an effective use of massive parallel mapping, partitioning and the necessary arithmetic 
operations. An example is offered for the simulation of control-structure interaction on 
a parallel computer and the impact of the approach presented herein for applications 
in other disciplines than aerospace industry is assessed. 


1. Introduction 

Active suppression of structural vibrations or active control of flexible structures has 
made considerable progress in recent years. As a result, it is now possible to actively 
suppress vibrations in mechanical systems emanating from machine foundations, in 
robotic manufacturing arms, truss-space structures and automobile suspension sys- 
tems. A common characteristic to these applications of active control theory has been 
its discrete actuators and discrete sensors, ranging from proof mass actuators and gyro 
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dampers to strain gages and accelerometers. Because most available discrete actuators 
are inertia force-oriented devices, actuation often triggers coupling between the actu- 
ator dynamics and structural transients. A practical consequence of such coupling is 
a limitation of achievable final residual vibration level if both the actuator and struc- 
ture possess insufficient passive damping level. It is noted that structures made of 
high stiffness composite materials have very low intrinsic damping, hence limiting the 
achievable residual vibration level for space maneuvering and space disturbance rejec- 
tion purposes. This has been a motivating factor for the development of distributed 
actuators and sensors which are often embedded as an integral part of the structure 
so that control force can be effectively maintained by strain actuation, thus alleviating 
the undesirable actuator dynamics associated with inertia-force actuation. 

Various activities that are being pursued by many investigators on the subject of 
adaptive structures may be categorized into three major thrusts: device developments, 
control laws synthesis and experimental demonstrations, and hard ware /soft ware im- 
plementation. The device developments effort has been the objective of many material 
scientists [1-3]. As the applications needs increase it is expected that functionally more 
reliable elect rostrictive and magnetostrictive elements will be available for use in active 
control/strain damping with improved product quality. 

The study of control laws synthesis and demonstration employing adaptive ele- 
ments has been one of the predominant activities in recent years. As scientists accu- 
mulate experience in the characterization of the coupling between the structure and 
the adaptive element, the applications will then be expanded from the current beam- 
like structures to the truss long beams, plates and shells. In order to effectively uti- 
lize as many adaptive elements as necessary for actively controlling the vibration of 
such large-scale structures in real-time operations, it will be imperative that the soft- 
ware/hardware components in the retd-time control loop must be able to process data 
fast enough so that control commands and the measurements can be carried without 
saturating and/or jamming the control system. 

With the advent of new technology in distributed actuators and sensors [4-9], it ap- 
pears that a combination of decentralized/distributed and hierarchical control strategies 
can be a viable alternative to conventional centralized control strategies. The real-time 
computer control of such systems as well as design of such control systems through 
iterating on simulations and hardware realizations thus will require the processing of a 
vast amount of data from and to the distributed actuators and sensors. A significant 
part of such data processing for the decentralized actuators and sensors is planned to be 
self- man aged, viz., there will be embedded microprocessors for each actuator and sensor 
pair or for each group of them. However, the necessary links between the decentralized 
control systems and the global control system as well as the necessary global control 
strategy will still require computational power far in excess of presently available real- 
time data processing capability. In addition, if one contemplates the performance of 
neural-network control or adaptive control for onboard real-time control of large-scale 
space structures, the computational need will dramatically increase beyond the current 
capability. As a case in point, even for the control of 20-bay truss beam vibrations by 
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three proof mass actuators and six sensors, NASA/Langley is relying on CRAY-XMP 
for adequate real-time data processing requirements. 

The objective of this paper is thus to present a computational framework by which 
one can bring the two emerging new technologies together, namely, the distributed actu- 
ators and sensors and the parallel computing capability, toward the real-time control of 
vibrations in large structural systems such as space stations, space cranes and in-space 
construction facilities. We will then discuss the potential for applying such a space 
technology to mitigate and/or minimize the earthquake damage of ground structures 
such as high-rise buildings, bridges and lifeline equipment. 


2. Models for Structures with Embedded Actuators and Sensors 

The coupling between the structural behavior and an adaptive electrostrictive element, 
whether it is embedded or surface-mounted, is primarily due to the following constitu- 
tive relation [3,10-12]: 


{:}- 

where e and v are the electrical displacement (charges/unit area) and the electric field 
(volt/unit area), <r and e sure the stress and strain, and s, g and c are the constitutive 
coefficient matrices, respectively. For magnetos! rictive elements, one needs to replace 
e and v by the magnetic field (H) and the magnetic induction (B), respectively, and 
the subsequent derivations will hold without any loss of generality. 



The coupled equations of motion for the structure and the adaptive elements 
can proceed by augmenting the standard procedure for the structure with the elec- 
tric transient equations plus the appropriate modification of the structural equilibrium 
equations that reflect the coupled constitutive equations (1). The resulting coupled 
structural-piezoelectric equations of motion take the following form [13-15]: 


( Structure: 


a) Mq + Dq + (K, + K a )q = f -i- Sa 

q(Q) = q 0 . q(0) = q 0 


Sensor Output: 

Actuator: 

Controller: 


6) y = H p q + H r q 4- H a a 

c) a+0a = B a u-S r |^ 

d) u + Gu = Ly 


(2) 




where 
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In the preceding equations, M is the mass matrix, D is the damping matrix, K, 
is the stiffness matrix due to structural strain-displacement relations and K„ is the 
stiffness matrix due to the strain actuation, f (t) is the applied force. S is the actuator 
projection matrix. H p , H r and H„ are the sensor calibration gain matrices, © is the 
actuator dynamic characteristics, B a is the gain matrix that translates the applied 
current /charge and voltage into the corresponding strain and strain rate where S is 
the transducer conversion gain, q is the generalized displacement vector and and the 
superscript dot denotes time differentiation, and u is the control law that consists of 
the applied current (or charge), Io, and voltage across the electrostrictive devices, Vo, 
G is the electric circuit characteristics, and L is the optimum direct feedback gain 
matrix. The case of dynamic compensations can be augmented to (2) by introducing 
an observer. But in subsequent discussions we limit ourselves to direct feedback cases 
only. 


It is noted that the control laws, unlike conventional control-structure interaction 
systems, are not directly fed back into the structural equations. Instead, the controller 
is simply a regulator controlling the electric charge, the voltage or the current. These 
regulated electric quantities are then fed into the piezoelectric sensors and actuators. 
Hence, it is the piezoelectric actuation that triggers feedback into the structures. 


3. Parallel Computations for the Dynamics of Adaptive Structures 

The earliest recorded computational results in mechanics were the parabolic trajectory 
calculations of a falling body by Galileo [16]. Since then, most scientific computations 
have been carried out by anthropomorphic algorithms, viz., step-by-step binary and/or 
decimal arithmetics. To set the stage properly for the present objective, parallel com- 
putations of the dynamic response of structures with distributed adaptive elements, we 
recall a passage by Kepler .to John Napier, the inventor of a logarithmic table: 

Newton was essentially dependent upon the results of Kepler's cal- 
culations, and these calculations might not have been completed but 
for the aid of that logarithms afforded. Without the logarithms, ...-, . 
the development of modern science might have been very different 

[ 17 ]. 

In terms of the present day data processing requirement, Napier’s logarithmic 
table in 1614 contained about 100 kilobytes, which was perhaps the most important 
computational aid to Kepler and Newton. Three and one-half centuries later we are 
witnessing gigabytes of tables being stored and retrieved at our disposal [18]. But 
these tables complement the weakness of the human mind and computational speed: 
long term memory and human arithmetic speed. In addition, for problems requiring a 
sequential nature of computations, i.e., ballistic trajectories which deal only with the 
position and velocity of a single shell or quasi-static equilibrium equations of a building 
structure, the computing activities do not interact with “time” and the computing 
efficiency affects only the humanpower efficiency for completing the computational 
task. - 
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There are many important scientific and engineering endeavors whose computa- 
tions must be fast enough for real-time delivery of the computed results. A classical 
example was Richardson’s lattice model for weather prediction by numerical process 
in 1922 [19]. The motivation for adopting such a lattice concept was due to the fact 
that the equation state at each lattice node takes on a different value set in time and 
an efficient way of interchanging and transmitting the nodal values at each time step 
was mandatory if the computations were to be carried out in real-time to predict the 
weather. Indeed, this was the dawn of the parallel computing era, even though the basic 
idea had to wait for its validity for 60 years. Today, many controls engineering activi- 
ties have been implemented by using computers so that their intended functions can be 
monitored and controlled in real-time. These include chemical processing, autopiloting 
and vibration control of simple structures. It is important to note that the computa- 
tional framework employed for such applications is based on sequential architecture. 
Hence, we believe that future improvements that can deal with large parameter models 
and large parameter controls must adopt a parallel computational framework. One 
such area is the dynamics and control of large structures with distributed/embedded 
adaptive elements. 

In order to carry out the necessary parallel computations, there are three distinct 
steps that must be addressed: discretizing the structure into appropriate partitions, 
mapping the physical partitions onto the processors, and step advancing of the equation 
states. These will be discussed below. 

3.1 Partitioning and Mapping of Adaptive Structures 

Ideally, if the sensor and actuator leads fall on the discrete nodes, no spatial interpola- 
tion would be necessary. However, such a situation is either difficult to realize or may 
prohibit the use of spatially convolving sensors [20] that are known to filter certain 
harmonic signals for minimizing phase lag in the feedback loop. Hence, we will assume 
that the sensor and actuator characteristics can be interpolated to the discrete nodes; 
in this way the partition boundaries can be chosen arbitrarily regardless of the physical 
locations of the sensor and actuator leads. In addition, this approach can lead to a 
natural embedding of the sensor and actuator characteristics into the finite element or 
boundary integral structural models. Once the partitioning is accomplished, the next 
step is to map the discrete partitions for adaptive elements onto the corresponding 
multiprocessors. 

Consider an adaptive structure that haw been modeled as a set of discrete elements 
as shown in Fig. 1. In a sequential computing environment, in order to advance the 
necessary computations for the present states, the arithmetic operations are carried 
out step-by-step for each node at a time. Hence, each nodal-state computations is 
performed in a manner similar to one courier delivering and picking up all the mails 
throughout the entire routes. In a parallel computational environment, in contrast, 
there can be as many couriers as necessary who comb through the routes concurrently 
in order to pick up and deliver all the mail at once. One of the most popular concepts 
in executing such tasks is the hypercube architecture (see Fig. 2) whose every node 
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Piezoelectric Elements 



Fig. 1 Discrete Model of Adoptive Structures 



Fig. 2 Hypercube Interconnection Network of a 
32-Processor 

(each node represents a processor) 
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is associated with a processor. Thus, to process the necessary computations for an 
adaptive structure with 19 partitions, one can assign the 19 adaptive elements to 19 
processors as shown in Fig. 3. The procedure for assigning the physical domain (ele- 
ments) to the parallel processors with minimal interprocessor communications is called 
mapping. 

Of several techniques available for the processor mapping of the computational 
domains [22], we will adopt a heuristic mapper developed by Farhat [23] since it can 
accommodate both the synchronous and asynchronous cases with robust and accept- 
able complexities. An application, of this mapping technique for modeling a bulkhead 
substructure for massively parallel computing is shown in Fig. 4. A similar mapping 
can be used for parallel computations of adaptive structures. 

3.2 Parallel Data Structure and Algorithms 

We will assume that each processor is assigned to carry out all the necessary computa- 
tions for at least one set of a sensor, an actuator, and a controller or a group of them. 
Therefore, the word partition does not necessarily imply a finite element: it can be a 
substructure, an element or even a sublayer within the composite layer that includes a 
sensor or an actuator. In carrying out the step-advancing in time, one may invoke an 
implicit or explicit direct time integration algorithm. When an implicit algorithm is 
employed, one needs to communicate not only the state variable vectors but also the 
associated matrices, i.e., the stiffness matrix, among the processors. Although we will 
show our results using implicit algorithms, we will, for illustrative purposes, restrict 
ourselves to an explicit direct time integration algorithm as it is intrinsically parallel 
and and the data structure aspects can be explained more succinctly via an explicit 
algorithm. It should, however, be mentioned that the choice of the solution algorithm 
can greatly influence the design and implementation of the necessary mapping and data 
structure. 

Consider the explicit integration of the equations of motion for the structure (2a) 
as recalled here: 

Mq -f f m4 = f -r f CO nt (3) 

where f,„< and f cont are the internal and applied control forces, respectively, given by 

fint = Dq + (K, + K a )q 

fcont = Sa 

The use of the central difference algorithm to integrate (3) leads to the following dif- 
ference equations in time 


X 


q" + * . q"“i + hM~ l (r + r eont - f? nt ) 
q n+ 1 = q" + hq n+ * 


(4) 
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(Physical Domain) 



(Processor Configuration) 


Fig. 3 Physical Domain and Its Mapping Onto 
Hypercube Processors 
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4. Implementation and Illustrative Example 


The mapping, partitioning and data structures above discussed have been imple- 
mented based on a shared-memory concurrent machine (Alliant FX/8) by modifying 
the software framework developed for finite element computations [26] and the control- 
structure interaction simulation and design software developed in (27, 28]. At present 
the following specialized systems of equations are implemented: 


where 


and 


Structure: a) Mq + Dq + Kq = f + Bu 4- Gw 

q(0) = q 0 , q(o) = q 0 

Sensor Output: b ) z = Hx +- m 

Estimator: c) x = Ax + Ef + Bu + L(z - Hx) 

x(0) = 0 

I Control Force: d) u = — Fx 


( 6 ) 


-{!}■ MI) 


H = [H, H.], L- , F=[F l F 2 


Li 


It is noted that in the above implemented equations, we have merged the actua- 
tor and the control law equations into one by neglecting the actuator and control law 
dynamics. Instead, we have introduced an estimator equation as we do not have all 
the measurements needed for complete feedback. In the above equations, B and ;3 
represent the input influence matrix for actuator locations whereas G and G represent 
the disturbance locations. The vector q is the generalized displacement, w is a distur- 
bance vector and the vector m is measurement noise. In Eq. (6b), z is the measured 
sensor output. The matrix is the matrix of displacement sensor locations and H„ 
is the matrix of velocity sensor locations. The state estimator in Eq. (6c) may or may 
not be model based. The superscript ' and " denote the estimated states and time 
differentiation respectively. The input command, u, is a function of the state estimator 
variables, q and q, and Fi and F 2 are control gains. The observer is governed by A, 
the state matrix representing the plant dynamics, and L, the filter gain matrix. 

The software thus implemented was used to test its applicability to solve the 
control-structure interaction design of a model Earth Pointing Satellite (EPS), shown 
in Fig. 6, which is a derivative of a geostationary platform proposed for the study of 
Earth Observation Sciences. Two flexible antennas are attached to a truss bus. Typ- 
ical missions involve pointing one antenna to earth, while tracking or scanning with 




Fig. G 


Earth Pointing Satellite Structure 
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Table 1„ EPS Vibration Frequencies (Hz.) 
Mode No. Frequency 


(1-6) 

0.000 

(7) 

0.242 

(8) 

0.406 

(9) 

0.565 

(10) 

0.656 

(11,12) 

0.888 

(13) 

1.438 

(14) 

1.536 

(15,16) 

1.776 

(17,18) 

3.026 

(19) 

3.513 

(20) 

3.531 


A small disturbance force was applied to the nominal EPS system in the form of a 
reboost maneuver. The force acted at the center of gravity in the Y-axis direction for 
0 seconds at a 10 N force level and from 0.1 to 0.2 seconds the force level was -10 N. 
The disturbance was removed after 0.2 seconds. Figure 7 shows the open-loop angular 
response about the X-axis of the 15 m antenna. A small amount of passive damping was 
assumed (D = 0.0002 K). The vibrational response produced more than 4.5 n radians 
of RMS pointing error due to this small reboost disturbance. Although many modes 
participate in the flexible body response, this particular reboost maneuver strongly 
excites modes near 4 Hz. The following paragraphs present an integrated control and 
structure design which seeks to lower the vibrational response of the EPS subject to 
some additional constraints. Figure 8 shows the closed-loop angular response about the 
X-axis of the 15 m antenna after design optimization. The pointing error is significantly 
reduced from that of the open-loop system shown . 


5. Future Work and Discussions 

The example problem analyzed in the previous section used a set of lumped actuators 
and localized sensors instead of distributed adaptive actuators and spatially integrated 
sensors. While such a model at best capture the adpative elements used by Anderson 
et al. [29], Matsunaga [30], and Takahara [31], it can not simulate on a large scale 
the distributed usage of piezoelectric actuators and sensors proposed by de Luis [32], 
Rogers et al. [33], and Burk and Hubbard [34]. Our immediate future work will 
concentrate on the implementation of distributed adpative elements and assess their 
practical applicability beyond the currently reported beam-like structural components. 
In this regard, we are exploring an adaptation of neural- network concepts [35] in the 
modeling and parallel computations of controlled structures with adaptive elements. 
Specifically, the limits of the applicability of distributed parameter modeling and control 
theory and discrete structures with discrete actuators and sensors, and their cross-over 
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Fig. 7 Open-Loop Response of EPS Structure 
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performance must be investigated. Design, modeling, simulation and testing criteria 
from such studies will provide greater insight into the eventual adoptions of adaptive 
structures as viable choice for future space systems design alternatives. 

The real-time simulation procedures presented herein may be applicable to the vi- 
bration control of lifeline equipment, and secondarily in minimizing the damage of build- 
ings during earthquakes. In this applications, the sensor measurements used herein can 
be directly applicable to the vibration and earthquake-causing forces on the structures. 
An idea that may prove to be crucial in this case is the use of earthquake-generated 
natural force as vibration minimization actuatori forces. In other words, instead of 
trying to mitigate the earthquake-generating forces, exploit the natural forces instantly 
to activate certain vibration minimizing devices! Research along this line may in the 
end lead to the design of actuators attachable to the columns and floors, if properly 
triggered during earthquakes, can minimize damages based on the natural forces. 
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